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Abstract 

An asymmetric double-well potential is considered, assuming that the minima of 
the wells are quadratic with a frequency lo and the difference of the minima is close 
to a multiple of Huj. A WKB wave function is constructed on both sides of the local 
maximum between the wells, by matching the WKB function to the exact wave 
functions near the classical turning points. The continuities of the wave function 
and its first derivative at the local maximum then give the energy-level splitting 
formula, which not only reproduces the instanton result for a symmetric potential, 
but also elucidates the appearance of resonances of tunneling in the asymmetric 
potential. 
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1 Introduction 



The quantum tunneling in a double-well potential appears in a variety of 
physical cases. Well-known examples include inversion of ammonia molecule 
[1], for which the double- well is symmetric, and proton tunneling in hydrogen 
bonds, for which the two wells could be unsymmetrical [2]. For a symmetric 
potential, instanton method developed in [3] , has been elaborated and applied 
to calculate energy splitting [4]. The tunneling splitting calculated from the 
instanton method exactly agrees with the WKB result when the quadratic 
connection formula is adopted, and it has been confirmed that the result is 
very accurate for large separation between the two wells [5,6,7]. 
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Quantum tunneling in asymmetric double-well potentials has also long been 
considered [8,9], and one of the intriguing properties in these cases is the 
appearance of resonances. For example, the wave function appropriate for the 
false vacuum of the potential 



[(x + a) 2 + (p 2 - a 2 )] for x < 



V D (x) = { V ' (1) 
^(x-P) 2 forx>0, 



(j3 > a > 0) can significantly tunnel to the side of x > only when (3 is tuned 
to satisfy the condition 



2 



((3 2 - a 2 ) « nhaj, n = 0, 1,2,3- ••. (2) 



Through various numerical calculations, it has been known that the appear- 
ance of these resonances is not limited to Vd(x), but is a general property of 
tunneling in asymmetric potentials [8]. In this paper, in order to elucidate the 
analytic structures of the resonances, we construct the WKB wave functions 
for a class of asymmetric double-well potentials. 

Specifically, we consider a smooth double- well potential V(x), assuming that 
V(x) has minima at x = b and at x = —a (a, b > 0), and a local maximum at 
x = 0. The minima are taken to be quadratic with a frequency u, V(b) = 0, and 
V(— a) = {n + e)hjj (See Fig. 1). As in the instanton method, we are interested 
in the large separation between the two wells, and consider the ground and 
low lying excited states of energy eigenvalue E satisfying V(0) ^> E > nftw. 
We also assume that the potential is still quadratic near the classical turning 
points between the wells. 

Around the minima, exact solutions to the Schrodinger equation are described 
by the parabolic cylinder functions [10]. On both sides of x = 0, a WKB wave 
function is constructed by matching the WKB function to the asymptotic 
forms of the exact solutions near the classical turning points. The continuities 
of the wave function and its first derivative at the local maximum then give the 
energy-level splitting formula. Though our method of requiring continuities is 
very different from the instanton or WKB method in [4], the splitting formula 
reduces to the known one in the symmetric case [4,11,12]. 

In the symmetric potential, for a given energy, an approximate solution to 
the Schrodinger equation localized in left (right) well implies, by the inversion 
symmetry, another solution localized in the right (left) well, and this fact has 
been conveniently used to evaluate energy splittings [13,11]. In this paper, we 
also show that tunneling in the asymmetric potential of e = can be explored 
by assuming the degenerate approximate solutions to the Schrodinger equation 



2 



V(x) 




Fig. 1. An asymmetric double-well potential V(x): V(b) = 0, V(—a) = (n + e)Huj, 
(n = 0, 1, 2, • • •). We assume that, for a given energy E, V(x) is quadratic in the 
regions of classical motions with the frequency u, and concentrate on the case of 
e< 1. 
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i/jr(x) and i/jl(x) which are localized in the right and left wells, respectively. 
By explicitly constructing iPr(x) and ^l(x), the splitting formula is found 
from these wave functions. Indeed it turns out that the WKB wave functions 
satisfying the continuities could be written as linear combinations of i/jr(x) and 
while the linear combinations lift the degeneracy to make the splitting. 
Therefore, a linear combination of the time- dependent WKB wave functions 
gives a system which shuttles back and forth between iPr(x) and ipi{x), to 
clearly elucidate the resonance structure of the tunneling in the asymmetric 
potential of e = 0. 

This paper is organized as follows. In Section 2, we construct the WKB wave 
function on both sides of the local maximum and, by requiring the continuities 
at the maximum, we evaluate the energy splitting. In Section 3, we find the 
appropriate localized wave functions 4>r(x) and iPl(%), to re-obtain the energy 
splitting formula. We also establish a time-dependent WKB wave function of 
a system shuttling back and forth. In Section 4, we give some concluding 
remarks. Finally in Appendix A we give exact solutions for the system of 
Vd(x) in the limit of large separation between the two wells. 



2 WKB method with continuity requirements 



As the results could be easily modified to incorporate the small non-zero e, we 
start with e = 0. 



2. 1 WKB wave function for x > 



For the wave function ipi(x) of the energy eigenvalue E = (i/+n+^)Huj around 
the right minimum, the Schrodinger equation is written as 

h 2 d 2 , mu 2 , t . 1., 

+ {x-byip! = hu){v + n + -)ip I . (3) 



2mdx 2T 2 v 2 
By introducing 



'Jl 

m/jj ' 



and zr = \[2{x — b)/lh a , we rewrite the equation as 



d 2 i) 



I , 1 ZRs , 

- + {v + n + -- -f 



dz\ ' ' " ' 2 



(4) 
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to obtain 



^{x) = C R D v+n {z R ) = C R D u+n ( ^y b) ), 



ho 



(5) 



where Cr is a constant and D v+n denotes the parabolic cylinder function [10]. 
Bearing in mind that we wish to construct a normalizable wave function, we 
choose the solution in (5) so that / 6 °° \ipi(x)\ 2 dx is finite if the expression of 
ipi(x) is valid for x > b. Asymptotic expansions of the parabolic cylinder 
function are well-known. For large and negative z (z <C —land z <C — \k\), we 
have 



D k (z)~e ~z k 



\ fc(fc-l) | 
2z 2 



T(-k) 



2vr u ■ z£ 
e e * z 



1 (fc + l)(fc + 2) 
2z 2 



(6) 



For x > 0, a classical turning point may be written as 
x = b u = b — \/2v + 2n + 1 Z/j . 



(7) 



In the classically forbidden region of 6j, — x ^> Z^, within the WKB approxi- 
mation, a solution to the Schrodinger equation is 



ip H (x) 



Ar p PM dv , B R _ [" pM 
p.Jo ft y A p Jo ft 



dy 



p(x) 



p(x) 



where p(x) is defined as 



p(x) = yj2m [V(x) - E], 



(9) 



and Ar, Br are constants. In the region of quadratic potential satisfying 
b — x ^> b — b u , by introducing 



p(y) 
h 



dy 



[ib-yf-{b- b v f] 

L ho 



21 1/2 



-dy, 



(10) 



we have [14,11] 



(b - xf 1 , 
<$>R(x) = - K -—^ + - A (2vAr2n + l) 



21 



ho 
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and thus 



A R ( 2^(b - x) 



v+n+i 



\Jmio{b — x) 



b-b„ 



exp 



B 



■s—(- 



6-6, 



^mu(b - x) V 2 Ve(6 - x) 



exp 



g(y) . _ (6 -xf 
n y ~ 211 



rp(y), . {b-xf 

dy + 



n/2 



(12) 



As we are interested in the limit of large separation between the two wells 
where the energy splitting is small, we introduce 5i with an integer I (= [u]) 
as 



6i = v -I, 



(13) 



so that 

\Si\ < 1. 



(14) 



In the region of the quadratic potential the wave function is also described 
by i>i{x). Making use of the asymptotic form in (6), for b — x ^> Iho, in the 
leading orders we obtain 



C 



R 



^ 'V2(x-b)\ l+n 



e ho 



+ 5,e 2l io 



2vr(/ + n)! 



>l,o 



V2l^(x-b) 



l+n+l 



(15) 



By matching the asymptotic form of ipi(x) onto that of i/ju(x) in this overlap 
region, we have 



.4 



R 



\l+n 



H(l + n)\g l+n _f»iEM dy 



I ho 



B R = (-l) l+n+1 C R 5 l 



27T 3 / 2 h(l + n)\ fh p(y) dv 

\ ; e Jo h , 

\ l' ho 9l+n 



(16) 
(17) 
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where 



WKB wave function for x < 



Since ^(x) = mw 2 (i + a) 2 /2 + n/tcj near the minimum of the left well, a 
classical turning point is given as 



x = —a. 



-a + \Zlv + 1 Z fto . 



(19) 



In the classically forbidden region of a v + x ^> 4, , a WKB solution may be 
written as 



/ ( \ A l -f° EM dy B L fo 

= e J* h y H ; eJ^ h y . 



p(:r) 



(20) 



where Al, Bl are constants. In the region of the quadratic potential satisfying 
a + x ^> a — a„, from the fact that 



$ L (x) 



ft 



/ 2 



(21) 



we obtain 



^m(^) ? 



a — a, 



5 L 



exp 



p(y) d _, (a + x) 2 



ft 



2/ 2 



2v^(a + x)V + ^ 



a — a. 



exp 



P(v)_ d _ (a + x) 2 



h 



o/2 



f22) 



With Zf, = v2(i + o)/lho, the Schrodinger equation around the left minimum 
is written as 



1 2 L 



+ iy + o - -r)^/v = o, 



2 4 
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and thus 

rl> IV (x) = C L D u {-z L ) = C L D u (- ^ + a) ), (23) 

<>ho 

with a constant Cl- The solution in (23) is chosen, so that \ipiv{x)\ 2 dx is 
finite if the expression of ?piy(x) is valid for x < —a. Around the left minimum 
of a + x ^> Iho, from (14) and the asymptotic form in (6), in the leading orders 
we have 



-lYCr 



6 h ° 



'V2(x + 



— 5i e 2l ho 



ko [V2l^(x + a) 

By matching ipiii( x ) to ipiv(x) in the overlap region, we obtain 



i+i 



(24) 



A L 
B L 



C L S l 



gJ— a; n " 

\ ko9i 



(25) 
(26) 



2.3 Continuity and energy splitting 



For a smooth potential, a wave function and its first derivative must be con- 
tinuous at x — 0, which gives the relations 



A) 



A 



Br = Br 



(27) 



There are three unknowns Cl, Cr, Si in these two equations, while an- 
other equation may come from the normalization of the wave function. From 
A L /B L = A R /B R , making use of Eqs. (16,17,25,26), we obtain 



5? 



9l 9l+n 



exp 



p(y) 
h 



dy 



-at 



{21 



which indicates that the splitting of the energy level A; is given by 



A; = yjg\ gi +n — exp 

7T 



p(y) 
h 



dy 



(29) 
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For a symmetric potential of n = 0, (29) exactly agrees with the result in 
[11,12]. 

Expression in (29) is not easy to use, because the integrand in the exponential 
is close to a singularity near the limits. By introducing 

b 

I a = J pm[V{y)-nhw]dy, I b = J y/2mV(y)dy, (30) 



and 



Vmw 2 1 , 

7o = / / \dy, 

'2 [V(y -a)- nhu] V 



o 

b 



pV(b-y) V 
the splitting is written for a general potential as 

A ,^^^(4^y +i/ 74^y + ° +i/2 . (32) 

^ir(l + n)\ l\ \ ho J V «fco / 



The expression in (32) can be conveniently used to find that our formula 
reduces to the known one in [12] for a symmetric potential. 



2.4 For a non-zero e 



The above formalism can be modified to include non-zero e, as far as 5i — e <C 1. 
In this case, without a change in ipi(x) and ^jj(x), the modifications of ipm(x) 
and ipiv{x) are obtained by replacing Si with 5j — e (or, 1/ with v — e). Due to 
the changes in Eqs. (24,25), the continuity requirements then give the relation 

W-) = (^) 2 . (33) 

If A^ denotes the energy splitting in the presence of e, (33) implies 

A^ = yjAf + (hue)*. (34) 
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3 An alternative method with localized wave functions 



For e = 0, if we assume two states of the normalized real wave functions i/jr(x) 
and with energy Eq as mentioned in Section 1, the Hamiltonian in this 

two- state subspace may be given by 



H 



E A/2 
A/2 E 



(35) 



where the small tunneling splitting A is written as 

h 2 d 2 



A = 2 



2m dx 2 



+ V(x) ipji(x)dx 



(36) 



The two eigenstates of the Hamiltonian are given by (4>r(x) ± ipi,(x))/ \[2. 
From the Schrodinger equations for these eigenstates, and the equations 



h 2 d 2 



2m dx 2 



+ V{x) 



i/ji(x) = E tpi(x) (i = R,L), 



and from the requirements 



(ijj R (x)) 2 dx w 1, 



{iPl(x)) dx ~ 1, / ip L (x)ip R (x)dx w 0, (37) 
o 



we find 



h 2 



A «-|Vl(0)^(0) -V*(0)Vl(0)| 



with ^(x) = dipi(x)/dx (i = R,L), which is a generalization of the method 
used for the symmetric case [11,13]. 



In the classically forbidden region, we may write 

iPr{x) = e J o h % Wl\x) = , e J o ft y , 



(39) 
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where Nr and Nl are constants. From these expressions of i/jr(x) and i/jl{x), 
by assuming that the validity condition for the WKB approximation 



d h 
dx p 



< 1 



(40) 



is satisfied at x = 0, we obtain 



A « 2- | N L N R 
m 



(41) 



For Eq = (I + n + |)ftu>, near the right-hand well, iPr(x) would be accurately 
described by the (I + n)th harmonic oscillator eigenf unction. This description 
holds well into the forbidden region, and, for (b — x)/lh 3> 1, we may write 



exp 



{x-b) 2 



2W 

no 



V2(x-b)\ 



l+n 



\JV* ko(l + n)\ \ l ho 



) 



(42) 



On the other hand, making use of (11), we can find the asymptotic expansion 
form of i[)r(x) of (39) in the overlap region. By matching the asymptotic form 
onto the expression in (42), we obtain 



N, 



R 



^y +n \/h gi+n 



2,71 I 



cxp 



h o 



p{y) 
h 



dy 



(43) 



Similarly, by matching the asymptotic form of i])l(x) onto that of the Zth 
excited harmonic oscillator state near the left-hand well, we have 



Nr 



2ix I 



■ cxp 



ho 



p(y) 
h 



dy 



(44) 



By plugging these explicit forms of Nr and Nl into (41), for E = (l+n+^)fiw, 
we confirm that A reduces to A/. 

If (27) is satisfied, after some algebra, we find that ipn(x) and ipni(x) of e = 
can be merged, in the classically forbidden region, into 



2ttA r I 



ho 



exp 



rp(y) 



h 



dy 



(45) 
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namely, iPwkb( x ) becomes 4>n(x) for x > 0, and ijjiii{x) for x < 0, where 
^Pwkb^) (^wkb^)) i s given when we choose 5i > (5; < 0) in (28). Thus 
this alternative method is in fact equivalent to the WKB method of requiring 
the continuities. We also note that iPwkb( x ) has a node in the classically 
forbidden region between the wells, while i^wkb{ x ) has no node in the same 
region. 

A (unnormalized) time-dependent WKB solution is given as 



\/2-KA R l ho 
2 exp 



-ioj{n + l + \)t + f^dy 
2 J n 
o 

A, 



-l) l+n cos (^j:t) Vfl(z) - i sin (^±t) M x 



(46) 



This last form shows clearly that the system shuttles back and forth between 
iJjr(x) and iPl(x) with the frequency Ai/h. 

In order to include e, let us consider a slight modification of the potential 
around the left well so that V{— a) changes from nftw to (n+e)ftw. While iPr{x) 
would still be an appropriate solution of the new system with energy .£/0) we; 
introduce ^£( x ) as an approximate solution with energy E + ehuj localized in 
the left well. If we confine our attention on the two state subspace described by 
iI>l(x) and i/jr(x), the Hamiltonian of the new system is analogous to that of a 
particle in a magnetic field [15] . We note that this analogy can be used to derive 
(34) within the approximation that ip e L (x) is the same with ip^x). If ip e (x,t) 
is a solution in this subspace with ip e (x, 0) = i/jr(x), the maximum of the 
probability of | j^ 00 {'ipl{x))*ip e {x, t)\ 2 during the time-evolution is A z 2 /(A z 2 +e 2 ), 
which indicates that the resonance peaks in tunneling have the Lorentzian 
shape. 



4 Concluding remarks 



Energy splitting formula has been obtained for the asymmetric double-well 
potential, by assuming that the potential is quadratic near the minima. As 
has been well known for the symmetric case, we expect that the splitting 
formula given here would be very accurate for the large separation between 
wells, which needs to be confirmed through numerical calculations. If we could 
add a linear term sx to the potential V(x) of e = with a controllable constant 
s, V(x) + sx has two minima at x = b — s/ muj 2 and at x — —a — s/ muj 2 . Since 
the difference of the minima is given as nfiw — s(a + b), in the light of numerical 
results [8] , the tunneling would be significant only if s is close to a multiple of 
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fku/(a + b). It would be of great interest to realize the asymmetric system with 
controllable constants. As a final remark, for n = and e = 0, we note that 
A 2Z /A = (2abe^ + ^/l 2 ho ) /(2l)\ which shows the quasi- Weierstrassian nature 
of the tunneling spectrum [12]. In this asymmetric case, thus, the tunneling 
behavior of an initially squeezed wave packet is erratic, and the trajectory of 
the expected position of the wave packet has a fractal structure. 



Appendix A 



For the potential Vd(x) in (1), WKB method may not be applicable, since 
the potential is not differentiable at x = 0. In this case, however, exact wave 
functions could be written in terms of the parabolic cylinder functions on 
both sides of x = 0. From the continuities of the wave function and its first 
derivative at x = 0, assuming m ^-{(3 2 — a 2 ) = (n + e)hw, we find that the 
eigenstate of an energy eigenvalue (v + n + e + \)hjj exists if the condition 



is satisfied. (47) can be solved in the limit of a, /3 > and e<l. Making 
use of the asymptotic expansion of (6), in this limit we obtain 

5? + (r(i2, - h) + e)5 l - R.U - erh x = 0, (48) 



where 



R, 



—— p " 1 ho 



(3 — a 
a + /3' 



2ir(l + n)\ 



\2l+l 
ho) -a 2 /l 2 h 

P. 1 ho _ 



2tt/! 



(49) 



(48) implies that the energy splitting is given as 

fiwy/AftU + e 2 + 2er(R l + U) + r 2 (^ - U) 2 . (50) 



If we formally use the formulas in Eqs. (30,31) by replacing V(x) with Vu(x), 
Ai coincides with the splitting of (50) in the symmetric case (n — 0, e = 0). 
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